Phase locking in quantum and classical oscillators: polariton condensates, lasers, and 

arrays of Josephson junctions 
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We connect three phenomena in which a coherent electromagnetic field could be generated: po- 
lariton condensation, phase-locking in arrays of underdamped Josephson junctions, and lasing. All 
these phenomena have been described using Dicke-type models of spins coupled to a single photon 
mode. These descriptions may be distinguished by whether the spins are quantum or classical, and 
whether they are strongly or weakly damped. 
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I. INTRODUCTION 

Phase-lockingpl of coupled oscillators is a well-known 
phenomenon in nonlinear dynamics. The generation of 
coherent radiation from Josephson junction arrays|2j is 
but one example for macroscopic oscillators. But phase- 
locking exists not only in classical systems but in quan- 
tum models. Here there are two basic paradigms for co- 
herence of microscopic oscillators: Bose-Einstein conden- 
sation (BEC), which is responsible for superfluidity and 
superconductivity, and lasing. 

Although BEC, lasing, and classical phase-locking all 
involve collective coherent behaviour, they are usually de- 
scribed in very different terms. With some exceptions ||, 
descriptions of BEC and lasing are given in quantum- 
mechanical language, which gives the impression that 
these phenomena derive from quantum mechanics. De- 
scriptions of Bose condensation sometimes go further, 
suggesting that the condensate itself is a quantum me- 
chanical object. Nevertheless, it is unclear how these 
phenomena differ from classical phase-locking. 

The aim of this paper is to clarify the relationships 
amongst BEC, lasing, and classical phase-locking, and 
hence the extent to which one can describe the coherence 
in Bose condensates and lasers as "quantum" or "clas- 
sical". To do this, we will consider Bose condensation 
of cavity polaritons, phase-locking in arrays of coupled 
Josephson junctions, and lasing. The simplest theoreti- 
cal models of these phenomena are in fact similar, and 
admit some controlled solutions, enabling us to compare 
them cleanly. 



be an unusual type of condensate. This is because low- 
energy polaritons are merely long- wavelength photons, 
which are not conserved particles. Thus the polaritons 
cannot condense in the ground state, making the bulk po- 
lariton condensate an intrinsically non-equilibrium phe- 
nomenon. However, a condensate of cavity polaritons is 
not necessarily a non-equilibrium phenomenon, because 
the lifetime of the low-energy cavity polaritons is finite. 
If this lifetime were long compared with the thermalisa- 
tion time, one could consider the thermal equilibrium of 
a population of polaritons, treated as conserved particles. 
This is the normal situation for Bose condensation. 

Eastham and Littlewood^, IToL ITT| have considered 
such quasi-equilibrium polariton condensation in a toy 
model. The model is appropriate to localised electronic 
excitations, such as excitons bound to impurities or lo- 
calised on quantum dots, dipole coupled to a single mode 
of a three-dimensional cavity. It has the Hamiltonian 



(1) 



ip' is the creation operator for a cavity photon, with en- 
ergy u> c . The dielectric is modelled as a set of N two- level 
systems, with the i th two-level system described by the 
spin-half operators Si . The eigenstates of S? correspond 
to the presence or absence of an excitation on site i. 

The model Q is the Dicke model^l °f quantum op- 
tics, which has been studied in many different regimes. 
To apply it to polariton condensation one should find its 
ground state or partition function, fixing the total num- 
ber of excitations 
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II. CAVITY POLARITON CONDENSATION 

A cavity polariton^, Q is the quantum of the electro- 
magnetic field in an optical cavity containing a dielec- 
tric. It is the confined version of the bulk polariton con- 
sidered many years ago by Hopfield|£|, which is formed 
from propagating photons coupled to electronic excita- 
tions such as excitons. Since polaritons are photons cou- 
pled to other excitations they are bosons, and therefore 
might be candidates for Bose condensation. 

While the idea of a Bose condensate of bulk polari- 
tons has been discussed for many years 0, 0], it would 
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which is conserved by the Hamiltonian Q). The con- 
straint on N cx can be dealt with by introducing a chemi- 
cal potential /i, so that one works with the unconstrained 
effective Hamiltonian i? c ff = H — fj,N ex . H c s is the same 
as H, except that the photon and exciton energies uj c and 
Ei are shifted by the chemical potential. 

One can write down the ground-state wavefunction of 
H e s by generalising the standard wavefunction for a Bose 
condensate. For bosons with creation operator b\ the 
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ground-state is the coherent state 



self-consistency condition on the synchronized states, 
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In general, polaritons are superpositions of an excitation 
of the cavity mode and an excitation of the dielectric. 
Thus the generalisation of J2J to describe polariton con- 
densation is 



exp(A^ 
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where A, Wi and <pi are variational parameters. Minimiz- 
ing (-ffeff) over these parameters gives an equation for A 



which is analogous to the BCS gap equation: 
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The polariton condensate © is a superposition of co- 
herent states of the dielectric and the electromagnetic 
field. It has a finite expectation value for the amplitude 
of the cavity field, (ip), and the electronic polarisation 
(S^). The cf>i are the phase differences between the elec- 
tronic polarisation and the cavity field. They are fixed 
by the dipole interaction term in which is responsible 
for the phase locking: it ensures that all the oscillators 
with a finite polarisation are mutually coherent, (pi = cp, 
when the energy is minimised. 

To make the connection to phase- locking more explicit, 
we note that the gap equation is a special case of the 
condition for the dynamics of the spins to synchronize 
at frequency [i. In a frame rotating at this frequency, 
the Heisenberg equations of motion corresponding to the 
Hamiltonian Q are 
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In an unsynchronized state, the sum on the right-hand 
side of Eq. [5] is of order y/N, so ip is of order 1 . For such 
a ip the spins are free to leading order in N due to the 
scaling of the coupling constant. The spin on site i sim- 
ply precesses around the z axis at its natural frequency 
{Ei — fi). In a synchronized state, part of the sum in Eq. 
El will be of order N, so ip will be of order y/~N. Such 
a field gives a finite contribution to the effective mag- 
netic field on each spin. The dynamics of (S~) in this 
effective field contains a static component. Substituting 
this static component into JSJ and setting ip = gives a 
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N ^ v /^-M) 2 +4 5 2 |A| 2 ' (8) 
This condition is a generalisation of (@J) , in which the unit 
numerator in the i th term of the sum becomes 2(S*')o. 
Sf is the component of spin i along its effective field, 
and ()o denotes an expectation value in the initial state. 
The particular choice of "initial conditions" that appears 
in Eq. 01 corresponds to a thermal (here T=0) occupation 
of the quasiparticle states, producing a solution of Eq. ^ 
of the lowest free energy. 

The self-consistent approach suggests that a mean-field 
theory in the amplitude of the cavity mode is exact as 
N — > oo. This can be formally demonstrated by con- 
structing the partition function as a path integral, which 
can be evaluated using saddle-point techniques. Physi- 
cally, the mean-field theory is exact because the cavity 
mode is coupled to many electronic states, and so should 
have relatively small fluctuations. A consequence is that 
condensation in the model Q is, in the limit of a large 
system, no more than the phase-locking of classical cou- 
pled oscillators. 



III. JOSEPHSON JUNCTION ARRAYS 

A different system which can be described by mod- 
els similar to £J is a Josephson junction array in a mi- 
crowave cavity. Phase-locking in that system was consid- 
ered in 1970 by Tillev[T^|. in a model describing junctions 
connected in series in a single-mode cavity. He consid- 
ered the fully synchronized states of the array, in which 
each junction oscillates at the same frequency. While 
there are many such states, differing in locking frequency, 
phase configuration and photon number, he showed that 
one particular state is selected by the driving current and 
cavity losses. 

To compare locking in Josephson junctions to po- 
lariton condensation, we use the angular momentum 
representation [l4) of a Josephson junction. Each junc- 
tion is represented by a spin S, whose magnitude is half 
the total number of condensed Cooper pairs in the two 
superconductors connected by the junction. S z is half 
the difference between the number of condensed pairs on 
either side of the junction, while S + and S~ transfer con- 
densed pairs across the junction. These operators obey 
the usual angular momentum commutation rules to the 
extent that Cooper pairs can be treated as structure- 
less bosons, which should be a good approximation in 
the weak-coupling limit A <C <-od- Note that S is usu- 
ally much larger than 1 — the junctions are themselves 
macroscopic objects. 



In the angular momentum representation, the canonical Hamiltonian for N Josephson junctions interacting with a 
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microwave resonance is 
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H = U! c ^1p + ^2 



Ki (Sf) 2 + (SJ 



ST)(il>-^)+3i(S+ + ST-) + Vi Sf 



(9) 



ji corresponds to the standard Josephson tunnelling, gi to the photon assisted tunnelling, Ki to the charging energy 
of the junction, and Vi to a voltage bias across the junction. To bring © as close as possible to the polariton 
condensation problem, we neglect the standard Josephson couplings ji and make the rotating-wave approximation. 
This gives the Hamiltonian 



(10) 



Since (|10|l conserves N ex , we may consider its behaviour at fixed iV ^[lfij. This problem should be exactly solvable in 
the limit N — > oo using a mean-field theory in the amplitude of the cavity mode. 

To construct the mean-field theory for H'jj = Hjj — /xiV ex , we should follow the same self-consistent approach that 
was used for polariton condensation. This is complicated, however, by the presence of the charging energy. Because 
of this term, the single-spin effective Hamiltonian of the mean-field theory is nonlinear. Instead of directly tackling 
this problem, we will follow the variational procedure used for polariton condensation. While this corresponds exactly 
to the mean-field theory for the polariton condensate, it will only be an approximation to the mean-field theory of 
the array. This is because © is the ground state of the effective Hamiltonian which occurs in the mean-field theory 
of polariton condensation, but not of that which occurs for the array. 

Looking for stationary points of (ifjj) in the variational state © gives equations for A and 0i — arctanu^: 



(w-/i)|A| = J^Bin(20i), 

i 

sin(2^)((y -v[) + 2K' l Scos{26 i )) = -2^|A| cos(2^). 



(11) 
(12) 



Here primes denote scaled variables, // = /iS etc., and the minimum energy solution has sin(20j) > 0. For some N ex 
and Ki, the charging energy term in l|12f) will be negligible. Eas.lllland ^]are then just the gap equation J3J derived 
for polariton condensation, with the replacements — > g' and \{Ej — fi) — > (y 1 — //). 



We have not investigated the consequences of the fac- 
tors of S and the charging energy in (|11H12|) . Never- 
theless, it seems that the mean-field theories for phase- 
locking in a Josephson array and for polariton condensa- 
tion are very similar. This is perhaps surprising, because 
we tend to think of Bose condensation as a truly quantum 
phenomenon occurring for microscopic oscillators, while 
Josephson junctions are macroscopic (S > 1), so that 
phase-locking is naturally thought of in terms of classi- 
cal coupled oscillators. However as we stressed in the 
last section, at zero temperature quantum mechanics is 
irrelevant to the mean-field theory of polariton conden- 
sation: the form of the gap equation (0} is the same for 
quantum spins and for classical angular momenta. The 
reason for this can be seen in the self-consistent dynam- 
ical approach, in which the problem is reduced to that 
of spins in a self-consistent field. For such a linear prob- 
lem, the commutativity or otherwise of the spin compo- 
nents is irrelevant. Notice however that there are at least 
two routes to the classical limit: in the Josephson ar- 
ray, the individual elements become classical as S — > oo, 
whereas for the polariton condensate we have S = 1/2, 
but N — > oo. In the latter system, only the coherent 
ground state can be treated as a classical object. 



IV. DECOHERENCE OF THE POLARITON 
CONDENSATE 

Some of the conspicuous differences between quantum 
and classical oscillators are due to the decoherence of 
quantum oscillators by their environment. Unlike a clas- 
sical oscillator, a quantum oscillator has states which do 
not have a well-defined phase. Furthermore, we expect 
that the environment will drive it towards such states. 
Thus we might expect that interactions with the environ- 
ment would have a significant effect on the phase-locking 
of quantum oscillators. In fact, we shall see that an infi- 
nite condensate is immune to weak decoherence phenom- 
ena, in the same way that a superconductor is immune to 
weak phase-breaking. But in the case of strong decoher- 
ence, and perhaps in the case of a finite system, we will 
find a connection to yet another example of macroscopic 
phase- locking — the laser. 

The laser and the polariton condensate are usually 
studied in separate contexts, and the connection between 
them is not made. This should be surprising, as both 
can be described by exactly the same Hamiltonian {Tjl. 
However, in a conventional laser |17| the only significant 
ordering is the coherence of the photons, while in the po- 
lariton condensate both the photons and the excitons are 
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coherent. This is because in a conventional laser the po- 
larisation of the gain medium, (S + ), is strongly damped 
by processes such as pumping, collisions, and interactions 
with phonons and impurities. The coherence in the pho- 
tons remains, because it can be generated by stimulated 
emmission even from an incoherent reservoir. 

The effects on the polariton condensate of different 
kinds of decoherence processes have recently been stud- 
ied by two of us [l8|, [l!J |2(| , using models related to . 
These models are obtained by rewriting each spin oper- 
ator in terms of a pair of fcrmions, with annihilation op- 
erators a, and bi. This is done by replacing with b\ai 
and Sf with \{b\bi — ajaj). With the local constraints 
b\bi + a| a, = 1 this would give an exact representation 
of the model Q). In our studies of decoherence, how- 
ever, we replace these local constraints with their global 
equivalent. The decoherence is modelled using baths of 
harmonic oscillators. Thus we consider the Hamiltonian 



H = H s + Hsb + H b . (13) 

The first term H$ is just the Hamiltonian Q written in 
terms of the fermionic operators, 




(14) 

Hb is a quadratic Hamiltonian describing the baths, and 
Hsb describes the coupling between the system and the 
baths. The most general form of Hsb is 

Hsb = Y,9 K k(^d k + dl^) 

k 

i,k 

+ if^ + a}a i )($+4) 

+ TfMh-ala^+ci)]. (15) 

The first term in (|15|) describes the decay of the cavity 
mode, the second term pumping of the two-level oscilla- 
tors, while the third term contains all the processes which 
destroy the electronic excitations, such as sponteneous 
emmission into modes other than the cavity mode. These 
baths could give rise to a flow of excitation through the 
system. The fourth and the fifth terms, however, describe 
all the dephasing processes which do not change the total 
number of excitations in the cavity, for example collisions 
and interactions with phonons and impurities. Processes 
described by the second, the third and the fourth terms 
in Ijl5|l have pair-breaking character, analogous to mag- 
netic impurities in superconductors, and correspond to 
potentials which vary rapidly in space or in time. 

In order to establish a crossover between an isolated 
condensate and a laser the decoherence processes must 
be included self-consistently. The widely used quantum 



Maxwell-Bloch (Langevin) equations with a constant de- 
cay rate for the polarisation are not correct when the 
coherent polarisation is large, i.e. for the polariton con- 
densate. In these equations the collective behaviour of 
the excitons is not taken into account when the lifetime 
for polarisation is derived. Instead, the lifetime for a 
single exciton is used in the equation for a collective po- 
larisation mode. 

To treat decoherence processes self-consistently, we 
use a procedure analogous to the Abrikosov-Gor'kov 
theory [2 lj of magnetic superconductors. In this theory, 
the baths which model decoherence are integrated out, 
introducing effective interactions between different two- 
level systems. These interactions are expressed as a self- 
energy in Dyson's equation, G,^ 1 = Gq ]j — Ey, of the 
form Sij = jGij. This form should be contrasted with 
the non self-consistent treatment, in which the decoher- 
ence appears as a constant lifetime in the Dyson equation. 

It turns out that the phase-locked polariton conden- 
sate is a robust phenomenon because at low decoherence 
strength it is protected by an energy gap proportional to 
the photon field amplitude. This gap becomes smaller 
and finally gets suppressed as the pair-breaking decoher- 
ence is increased. At low excitation densities this leads 
to the suppression of all the coherent fields while at high 
densities it leads instead to the conventional character- 
istic of a semiconductor "laser" in which the coherence 
is almost entirely in the photon field and there is no gap 
in the excitation spectrum. The laser regime of a polari- 
ton system emerges in a way that demonstrates its close 
analogy to a gapless superconductor. 

Although the coherent polarisation in a conventional 
laser is strongly damped, it must be finite for the laser 
to operate. Thus the transition between a condensate 
and a laser is smooth. There is no formal distinction 
between the two based on the broken symmetry of the 
ground state. There may be useful practical distinctions 
however, such as the presence or absence of a gap in 
the excitation spectrum. We do not yet know if the dy- 
namics of the order parameter, and hence the linewidth, 
differs. Note that in the Abrikosov-Gor'kov theory the 
pair-breaking does not produce fluctuations in the order 
parameter. The theory presented here may share this 
feature, whereas real lasers have a finite linewidth. 



V. FINITE-SIZE FLUCTUATIONS 

The question of classical or quantum behaviour can 
never arise for the dynamics of the order parameter in 
an (infinite) system with a broken symmetry. One sim- 
ply has a macroscopic equation of motion for the order 
parameter interacting with an external classical field — 
with a familiar example being the Josephson equation 
for a weak link, and a less familiar one the classical dy- 
namics of the (averaged) mean field equations in section 
II. However, when such a system becomes finite (though 
still large) in extent, we can ask whether the dynam- 
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ics are now best descibed by a Schrodinger equation or a 
Langevin equation: the order parameter will "diffuse" , at 
short times following quantum mechanics, and at longer 
times dictated by Brownian motion. Of course, we are 
now concerned only with the low energy degrees of free- 
dom — those near the frequency \x in the rotating frame 
— and certainly well away from the quasiparticle excita- 
tions above the gap. 

The procedure to be followed is clear, at least in prin- 
ciple, though it has not been fully completed for the 
model of a polariton condensate. As we mentioned briefly 
above, the variational equations correspond to the sad- 
dle point of a quantum mechanical action, which is exact 
as N — > oo. Fluctuations at finite N are described by 
a new effective action, with degrees of freedom that are 
then coupled to baths exactly as in EqEl However, in 
contrast to the results of the last section, phase-breaking 
perturbations are expected to be always relevant how- 
ever weak. We note in passing that this methodology is 
different from the conventional procedure to begin with 
a classical action that is then re-quantised. Whether or 
not it yields any distinct difference is not known. 

For a finite system, the broken symmetry ground state 
will not be stable, and we will observe fluctuations; this 
is a familiar point of view in classical laser theory, where 
Haken has emphasised how the mean field theory corre- 
sponds to a second-order phase transition, and the fluc- 
tuations in a real system arise because the number of 
photons is not infinite. However, with the starting point 
of an effective quantum mechanical action one no longer 
presupposes a classical limit. With the new action one 
can compute, for example, correlation functions of the 
photon field, which are directly measurable. At the mean 
field level, the photon field is a classical electromagnetic 
field, so the distinction between quantum and classical 
statistics will only appear at this stage. 



VI. DISCUSSION 

The distinctions between coupled oscillators, BEC of 
polaritons, and strongly damped lasers is an important 
one in the context of recent experiments on semiconduc- 
tor microcavities. We will not give a detailed review of 
the field here, except to point out that one recent ex- 
periment has shown evidence for coherence of the 
photon field in the nonlinear but incoherently pumped 
microcavity. However, if we define BEC of polaritons 
to be restricted to systems where the excitonic degrees 
of freedom have strong coherence, observation of photon 
coherence is by itself not decisive evidence of polariton 
BEC. One characteristic of this regime would be a gap in 
the excitation spectrum, which was not apparently ob- 
served. 

The difficulties of making these distinctions had al- 
ready arisen in the Josephson junction array problem — 
here embedded in microwave cavities. As we discussed 
in section IIIII the Hamiltonian is mathematically sim- 



ilar to that of the polariton condensate. Recently [23|, 
Barbara et al. reported thresholds in the ac output 
power of such an array as the driving power was in- 
creased. They interpreted their results as analogous to 
lasing, with gain due to stimulated Josephson tunnelling. 
Subsequently, Stroud et al. showed that many of the ex- 
perimental observations could be reproduced by classical 
treatments \24. l25t \2dL |27 | of models similar to 

The phenomena of lasing and BEC are well-known to 
be closely connected. At a microscopic level, both are 
consequences of the quantum mechanics of indistinguish- 
able bosons, and more specifically of stimulated scat- 
tering. Lasing is described dynamically, and the role 
of stimulated scattering is explicit. Bosc condensation 
is described thermodynamically, with the role of stimu- 
lated scattering hidden in the Bose-Einstein distribution. 
On a more sophisticated level, Haken |2{| has shown that 
the mean-field theory of the laser is analogous to that of 
a second-order phase transition, while Oraevskii[29| has 
discussed the dynamics of a superconductor in terms of 
stimulated scattering. 

What confuses the issue about BEC is that the conven- 
tional textbook picture presents BEC as a consequence of 
statistical physics of weakly or non- interacting bosons, 
which obscures the central point that BEC in a macro- 
scopic system is a phase transition like any other. So 
in a very large system one will not expect to find quan- 
tum mechanics operating at the level of the macroscopic 
order parameter, even if the microscopic theory of this 
phase transition requires quantum physics. In the JJ ar- 
ray, it is usual to imagine that the individual elements 
are macroscopic ("decohered") from their environment, 
but one can see from the above that is is not necessary 
to assume this in order to develop the correct classical 
theory of the phase locked coupled array. For our sim- 
ple model of polariton BEC one can see again that in the 
large system limit with macroscopic occupancy, the quan- 
tum mechanical ground state corresponds to the dynam- 
ics of classical phase-locked oscillators. But here there 
is a possibility to decohere the individual elements (spin 
half dipoles) from each other (by coupling to external 
baths) and restore conventional laser theory with a co- 
herent photon field, supported by incoherent electronic 
polarisation. But still, in the limit of macroscopic occu- 
pation, the coherent photon field is essentially a classical 
one. 

Most interesting would be the behaviour of large, but 
not infinite, systems. Here the order parameter fluctuates 
generically due to finite size effects. These fluctuations 
can arise due to environmental interactions, which will 
give rise to classical diffusion, as in the standard theory of 
the laser linewidth near threshold. They might also arise 
because the order parameter tunnels between equivalent 
states, as has been achieved in systems of small Joseph- 
son junctions[33,|^|. The competition between quantum 
mechanical tunneling and environmental dephasing is of 
course at the heart of current attempts to create quan- 
tum coherent devices — and if excitonic or polaritonic 
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BEC were observed, this would provide another possible 
fundamental system upon which to base such work. 
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